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ABSTRACT 


Theoretical aspects of a new capability developed and implemented in NASTRAN 
Level 17.7 to analyze forced vibration of a cyclic structure rotating about its 
axis of symmetry are presented. Fans, propellers, and bladed shrouded discs of 
turbomachines are some examples of such structures. The capability includes the 
effects of Coriolis and centripetal accelerations on the rotating structure which 
can be loaded with: 

1) directly applied loads moving with the structure and 

2) inertial loads due to the translational acceleration of the axis of 
rotation ('base' acceleration). 

Steady-state sinusoidal or general periodic loads are specified to 
represent: 

1) the physical loads on various segments of the complete structure, or 

2^ the circumferential harmonic components of the loads in (1). 

The sinusoidal loads are specified as functions of frequency and the general 
periodic loads are specified as functions of time. The translational acceleration 
of the axis of rotation is specified as a function of frequency in an inertial 
coordinate system. 

The cyclic symmetry feature of the rotating structure is used in deriving 
and solving the equations of forced motion. Consequently, only one of the cyclic 
sectors is modelled and analyzed using finite elements, yielding substantial 
savings in the analysis cost. Results, however, are obtained for the entire 
structure. A tuned twel ve-bladed disc example is used to demonstrate the various 
features of the capability. 


1 . INTRODUCTION 


mtrr 


') 


Under the sponsorship of NASA's Lewi". Research Center, a series of new 
capabilities has been developed and added to the general -purpose finite-element 
structural-analysis program NASTRAN [1-7]. A variety of problems including 
static aerothermoelastic and dynamic aeroelastic analyses of tuned cyclic 
structures, and modal analysis of mistuned cyclic structures, such as bladed 
discs of turbomachines, and advanced turbopropel lers have been addressed. 

This paper presents the theoretical aspects of one of these NASTRAN capabilities 
[4,5]. 

Figure 1 illustrates the problem by considering a 12-bladed disc as an 
example. The bladed disc consists of twelve 30° segments--identical in their 
geometric, material and constraint properties. The disc rotates about its axis 
of symmetry at a constant angular velocity. The axis of rotation itself is 
permitted to oscillate linearly in any given inertial reference. In addition, 
the bladed disc is allowed to be loaded with steady sinusoidal or general periodic 
loads moving with the structure. Under these conditions, it is desired to deter- 
mine the dynamic response {displacement, acct ation, stress, etc.) of the 
bladed disc. 
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2. EQUATIONS OF MOTION 


The cyclic symmetry feature of the rotating structure Is utilized in deriving 
and solving the equations of forced motion. Consequently, only one of the cyclic 
sectors is modelled and analyzed using finite elements, yielding substantial 
savings in the analysis cost. Results, however, are obtained for the entire 
structure. The Coriolis and centripetal acceleration terms have been included. 

For clarity of derivation, the equations of motion are first derived for an arbi- 
trary grid point of the cyclic sector finite element model, and then extended 
for the complete model. 

COORDINATE SYSTEMS 

These are shown in Figure 1. 0-XYZ is an inertial coordinate system. 
0-XgYgZg is a body-fixed coordinate system such that QXg coincides with the axis 
of rotation of the structure and is always parallel to OX. For a NASTRAN finite 
element model of the bladed disc, 0-XgYgZg also represents the Basic coordinate 
system. A-xyz is a body-fixed global coordinate system in which the displace- 
ments of any grid point P are desired. The unit vectors associated with these 
coordinate systems are also shown in Figure 1 . 

DEGREES OF FREEDOM 

W 

The rotating structure is permitted four rigid body motions including three 
translations (along f)X, ®Y and CZ) and one rotation at a constant angular velocity 
£2 about its axis of rotation OXg. 
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Each grid point of the structure is permitted six degrees of freedom. 

The displacement at any grid point in any sector can be expressed in any body- 
fixed coordinate system as a combination of: 

1) the steady state displacement due to the steady rotation of, and the 
steady state loads applied to, the structure, and 

2) the vibratory displacement (superposed on the steady displacement) due 
to the vibratory excitation provided by the directly applied loads and the 
inertial loads due to the acceleration of the axis of rotation ('base' 
acceleration). 

The purpose of. the present development is to determine the vibratory 
response. 

LA GRANGE FORMULATION 

Referring to Figure 1 , the complete tuned structure consists of N identi- 
cal cyclic sectors. If u represents all the vibratory degrees of freedom of the 
complete structure, the equations of motion can be derived via the Lagrange 
formulation, 


HxlU4,3D_9W 

dt 1 aft ' " 9u 3u aa ~ au * 


where T and U represent the kinetic and strain energies, respectively, of the 
complet e structure; D is the Rayleigh's dissipation function representing the 
energy lost in the system due to resisting forces proportional to velocities 8 
(e.g. viscous damping forces); and 6W represents the virtual work done on the 
structure by the external forces through virtual displacements 6u. 


The complete set of degrees of freedom u can be subdivided into N subsets, 
each containing u n degrees of freedom for each of the N cyclic sectors. Since 
any given cyclic sector is 'connected' to adjacent cyclic sectors only on its 
two sides, u n satisfies the intersector boundary compatibility condition 


u n s u n+l 
u side 2 u side 1 


n - 1 , 2, 


N. 


( 2 ) 


Equations (1), therefore, can be written as N sets of equations coupled 
only as given by equations (2): 


d / 3T n N 9T n , 9U n . 3D n . 3W n 

dt ( ^ ) “W r+ a^ + '§5 fr “'^ Tr » 


N . 


n = 1 , 2, . . . , 


( 3 ) 


For clarity of presentation, without loss of generality, equations (3) are 
first applied to obtain the equations of motion of an arbitrary grid point in any 
cyclic sector by considering its three translational degrees of fre?d;>m. Inclu- 
sion of the three rotational degrees of freedom at the arbitrary grid point, and 

extension to include the remaining grid points in the cyclic sector are con- 

* 

sidered subsequently. 


EQUATIONS OF FORCED MOTION 


With reference to Figure 1, point P is an arbitrary grid point of the nth 
cyclic sector with a mass of 'm' units lumped from the adjacent finite elements. 
Substitution of the expressions for T, U, D and 6W in the Lagrange equations (3) 
results in the following equations of forced motion of point P expressed in the 
displacement (global) coordinate system A-xyz: 


[M] fti} + 


[B] + 


2f2[B-j] 


(u) + 


[K] - srCMj] 


(u) = {?} 


[m 2 3 {R 0 > 

(4) 


T he terms appearing in equations (4) are given in Appendix A. 


Equations (4) describe the translatory motion of an arbitrary point P in 
an arbitrary sector n of the rotating cyclic structure subjected to a directly 
applied vibratory load {P} and base acceleration { R Q } . 

These equations can be extended to include the three rotational degrees of 
freedom at point P by noting that: 


1) in a lumped mass model, only the translational degrees of freedom at 
any grid point contribute to the kinetic energy of the structure, and 

2) the coupling between various degrees of freedom may exist only via 
the stiffness matrix. (Instances where the damping matrix is defined 
proportional to the stiffness matrix also may result in coupled equations 
of motion. ) 


by 


Accordingly, the matrices derived from kinetic energy considerations, 
[M], [B-j], [Mj] and [M 2 ] of equation (4 ) can be expanded as typified 


M 0 
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where [M tt ] Is the 3x3 (translational) mass matrix of equation (4 ). With sub- 
scripts t and r representing the translational and rotational degrees of freedom 
at point P, the stiffness and damping matrices may be expanded as 


CK3 


and 


[B] 


E^tt^ I ^tr^ 

[K rt ] 


t B rt-^ 


CK rr 3 


t_ B t£ 

CB rr ] 


( 6 ) 

(7) 


By similar reasoning, the equations of forced vibratory motion of all the 
cyclic sectors of the total structure can be written as 


[M n Ku' n } -v [;B n ] + {u n > + [tK n ] - f) 2 [M"j]{u n ) = {P n } - [M^J{R }, 


n * 1, 2 N. (8) 

The intersegment boundary compatibility is specified by equation (2). 
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3. SOLUTION OF EQUATIONS OF MOTION 


The method of solution of the equations of forced motion (equations 8 
and 2), is based upon the form In which the excitation of the rotating struc- 
ture is specified. Because of its eventual implementation in trie naSTRAN 
general purpose finite element structural analysis program, the following 
solution procedure is generally similar to the theoretical presentation of 
cyclic symmetry given in the NASTRAN Theoretical Manual [8], 


METHOD OF SOLUTION 

The method of solution of the equations of motion consists of four principal 
steps: 

1) Transformation of applied loads to frequency-dependent circumferential 
harmonic components. 

2) Application of circumferential harmonic-dependent inter-segment com- 
patibility constraints. 

3) Solution of frequency-dependent circumferential harmonjc components of 
displacements. 

4) Recovery of frequency-dependent response (displacements, stresses, loads, 
■ etc.) in various segments of the total structure. 
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An overall flowchart outlining the solution algorithm is shown in Figure 2. 
Provision to Include the differential stiffness due to the steady loads is also 
shown. 

1 . Transformation of Applied Loads 

♦ 

The transformation to frequency-dependent circumferential harmonic components 
depends on the form in which the excitation is specified by the user. 

Details of the- five loading conditions considered are as follows: 


Directly applied loads (segment-dependent and periodic in time) 

If P n represents a general periodic load on sector n specified as a function 
of time at M equally spaced instances of time per period (Figure 3 ), the load 
at m^ time instant can be written as 


m 


P n = P n + l 


*■1 


-£c -£s 

P n cos(m = TAb) + P n sin(n>T£b) 


•M/2 


+ (-1 ) m ”‘ l p n 


(9) 


m - 1, 2, • . . , M i 

where b * 2 tt/M, £. * (M-l)/2 for odd M, £, * (M-2)/2 for even M. The last term 

L L -"P," 

in equation (9) exists only when M Is even. The coefficients pn 

("•£" = 0; £c, £s , £=1 , 2, ..., £ l ; M/2) in equation (9) are Independent of time, 

and are defined by the relations 

-o M m 

P n = l l P n , (* ■ 0) Part of (10) 

M m=l 
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1 r/'f 

b c . _.a * J 


-Ac M m 

P n - | I P n cos(ffhTAb), 


o w§:m» : 

OF POOR OlwU iV 


- As 


ID*1 

M 


m 


y ( A*1 » 2, nil Ai ) 


pf1 * I pn sin(iivTAb), and 
m*l 


> (10 Contd, ) 


-M/2 


M 


m 


P * i I (-l) m 1 P n (M even only) (A-M/2). 

” jtf! 

-"A" 

Each of the coefficient vectors P n on the left hand sides of equations 
(10) can further be expanded in a circumferential (truncated) Fourier series 

. h o *• -’’p," r >£' 1 Mpii t -"a 11 

■ P° + T Ip * 6 cos (rPTka) + P« slnfiCTka) 


i. 


1 l- 


where n * 1 , 2, . . , , N , 

"A" * Oj Ac, As, A B 1, 2, . Jt L ; M/2 
a * 2 tt/N 

k, ^ (N-D/2 for N odd 


+ (- i ) n ‘¥ /2 , 

ni) 


( 12 ) 


k L =(N-2)/2 for N even. 

The last term in equation (11) exists only when N is even. The Fourier 

-V 

—"k" 

coefficients P ("k" =0; kc, ks, k = 1, 2, .... k L ; N/2) in equation (11) 
do not vary from sector to sector, and are defined by 


" 0 « 


-"A 


N 


-''A" 


T° *Jr I P" 


N 


n=l 


(k * 0) 


-"A" 


-"A" 


T* c -£ I P n cos (n^Tka) 


n=l 


-“A" 

■prks 


c" 

pN/2 


N 


fl P 

N n*l 


N 


-"A” 

n sin (n-lka), and 


(k *= 1, 2, ...» k L ) 


> ( 13 ) 


i'l .(-l)"-’ p" 


n=l 


(N even only) (k « N/2) 


ORIGINAL PASS JSf 
OF POOR QUALrrY 

The terms F ^ ("A" ■ 0; Ac» as, A « 1, 2, A L ; M/2 and *'k" * 0; kc, ks, 
k • 1, 2, . k L ; N/2) are the transformed frequency-dependent circumferential 
harmonic components of the directly applied loads (m * 1* 2, . M and n * 1 , 
2 N). 


Directly applied loads (Circumferential harmonic-dependent and periodic in 
time) . 


Such loads can be represented as 
m -o A, 


r"k M 


-0 

n"k" 


+ L L 


-Ac 

_JI L II 

P K 


-As 


cos(nPTAb) + p' k "sin (iTAb) 


-M/2 
(-1 )m _ lp"k" 


(14) 


where m * 1, 2, .... M represent the time instances at which harmonic components 

"k" ■ 0; kc, ks, k * 1, 2, ...» k, ; N/2 of directly applied loads are specified. 

-"A" L 

•The coefficients ^"k" on the right hand side of equation (14) are obtained 
using equations (10) with sector number n replaced by harmonic number "k". 

Directly applied lo^ds (frequency-and segment-dependent) 

This type of loads can be represented as 


-ll£ll _ll£H -1 „llpll 

P* c cos(rTTka) + F* 5 sin(ivfka) + (-l) 11 ' 1 P^ /2 , (15) 

where "A" (=1, 2, ..., F) now represents the frequencies at which excitation is 

specified. The transformed frequency-dependent circumferential harmonic components 
-"A” 

— "k" 

P ("k" = 0; kc, ks, k = 1, 2, .... k^; N/2) are obtained using equations (13) 
with "A" as defined above. 

Directly applied loads (frequency-and circumferential harmonic-dependent) 

These loads are the transformed frequency-dependent circumferential harmonic 
components ("k" * 0; kc, ks, k * 1 , 2, . . . , k L ; N/2) with "A" (=1, 2, ..., F) 

representing the various frequencies at which the directly applied loads are 
specified. 

Base acceleration (frequency- and circumferential harmonic-dependent) 

i. 

In Appendix B, it is shown that the components of the translational base 
acceleration contribute to inertial loads on the rotating structure in the follow- 
ing manner: 


-"A" -"A" k. 

P n * F° + l 
k R l 


-n- 


1 . 


-"A" 

it b n 

Axial component contributes to F * where "k" * 0, and "Jl" represents the 
specified excitation frequencies, 

— "k" 

2. Lateral components contribute to P where "k" * 1c and Is, and re- 
presents the ef fective excitation frequencies which are shifted from the 
specified frequencies by i ft, the rotational frequency. 


2. Application of Inter-Segment Compatibility Constraints 

As shown in Section 4.5.1 of Reference 8, equations (2) 
are used to derive the compatibility conditions relating the circumferential 
harmonic component degrees of freedom on the two sides of a rotationally cyclic 
sector: 

side 2 side 1 

(k « o) 


—0 

u 2 


—0 

U 1 


uf * u^ c cos(ka) + u^ s sin(ka) 


U 2 S = -u^ c sin(ka) + u^ 5 cos(ka) 




> (k * 1, 2 k L ) > (16) 


and 


u ^ 2 


-uj ' 2 


(k = N/2) 


In order to apply these constraint relationships for any given harmonic k, 

_!/ 

an independent set tr consisting of the circumferential harmonic component (cosine 
and sine) degrees of freedom from the interior and side 1 of the cyclic sector is 

i/ 

defined. Ir is selected from the 'analysis' set degrees of freedom (i.e., the 
degrees of freedom retained after the application of constraints and any other 
reduction procedures), and is defined as 


u* c = G ck (k) i^, and 

u* s * G $k (k) u* 


> (17) 


u* c and u* s each contain all (and only) the 'analysis' set degrees of freedom 
from the interior and both sides of the cyclic sector. Equations (16) are used 
to define some of the elements of the transformation matrices G ck and G sk « For 
k * o and N/2, the matrix G $k is null. 


ORIGINAL PA 1 :) I I J 
OF. POOR QUAUYV 


3 , Solution of Frequency-Dependent Harmonic Displacements 

1/ 

For a given harmonic k, the introduction of iT in the equations of motion, 
equations (8), results in the transformed equations of motion 


where 



+ 1*5* 

i/ y t/ 

+ « r 

« 


! ^ * 

r T 

6 ck 

M n 

G ck + G sk Mn 

G sk • 


I* « 

r T 

G ck 

B n 

G ck + G Ik B " 

G sk > 


K* * 

r T 

G ck 

K n 

G ck + G Ik *" 

G sk * 

and 

■pK E 

r T 

G ck 

pkc + G T k 

1 


As discussed 

earlier, P* 0 

and 

pks 


(18) 


(19) 


frequency-dependent circumferential harmonic components of the directly applied and 
base acceleration loads. 


At any Citation frequency o/, let 
= P K e^ t and accordingly, 


-K =K iw"t 
IP = u e 


( 20 ) 


where P K and u^ are complex quantities. Equation (18) can be rewritten as 

[-u.' 2 # + iu,'B* + K K ]U K = P K . (21) 


i 


( 22 ) 


The excitation frequency w' is given by 

w" * w for all directly applied and axial base acceleration 
loads, and 

■ w±n for lateral base acceleration loads. 

* 

Equation (21) is solved for u for all excitation frequencies and all harmonics 
as specified by the user. The cosine and sine harmonic components of displacements 
are recovered using equations (17). 

4. Recovery of Frequency- Dependent Displacements in Various Segments 

This step is carried out only when the applied loads are specified on the 
various segments of the complete structure. 


For loads specified as functions of time, equation (11) is used to obtain the 
displacements In various segments with « 0; Ac, As, i * 1 , 2, .... 

u max 

The circumferential harmonic k Is varied from k^ n to k max . The user specifies 

_ .. 3 k M j _ an d k . . . y • 

max mm max 

For loads specified as functions of frequency, equation (15) is used to 

obtain the displacements u n in various segments with representing the 

excitation frequencies. The circumferential harmonic is varied from user speci- 

fied k . to k , 
mm max 

The theoretical development discussed above has been implemented as a new 
capability in NASTRAN Level 17.7. All aspects of the capability, including 
DMAP changes, are extensively documented in References 4 and 5. 


4. EXAMPLES 


Five inter-related examples are presented to illustrate the theoretical 
development of the previous sections. The new capability added to NASTRAN to 
conduct forced vibration analysis of rotating cyclic structures [5] has 
been used to conduct these examples. A 12-bladed disc is used for 
demonstration. 


Example 1 is conducted on a finite element model of the complete structure 
(Figure 4). Examples 2 through 5 use a finite element model of one rotationally 
cyclic sector (Figure 5). Results of example 1 are used to verify some of the 
results obtained in the remaining examples. Table l summarizes the principal 
features of these examples. 


Steady-state frequency-dependent (sinusoidal) or time-dependent (periodic) 
loads are applied to selected grid point degrees of freedom. The specified loads 
can represent either the physical loads on various segments or their circumfer- 
ential harmonic components. For illustration purposes only, the frequency band 
of excitation, 1700-1920 Hz, due to directly applied loads and base acceleration 
is selected to include the second bending mode of the disc for a circumferential 
harmonic index k * 2. The 'blaae-to-blade' distribution of the directly applied 
loads also corresponds to k = 2. Table 2 lists the first few natural frequencies 
of the bladed disc for k = 0, 1 and 2. Modes for k = 2 are shown in Figure 7. 


GENERAL INPUT 

1. Parameters: 


Diameter at blade tip s 

Diameter at blade root s 

Shaft diameter = 

Disc thickness s 

Blade thickness * 

Young's modulus = 

Poisson's ratio s 

Material density * 

Uniform structural damping 


19.4 in. 

14.2 in. 

4.0 in. 

0.25 in. 

0.125 in. 

30.0 x 10 6 lbf/in 2 
0.3 

7.4 x 10“ 4 lbf-sec 2 /in 4 

(g) * 0.02 
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2. Constraints: 


A cylindrical coordinate system with origin at the Basic origin, Z axis 
along Basic X axis, and 9 measured from C;e Basic Y axis is chosen as the global 
(displacement) coordinate system for all grid points. 

All constraints are applied in body-fixed global coordinate system(s). 

All grid points on the shaft diameter are completely fixed. Rotational degrees 

of freedom 0 7 at remaining grid points are constrained to zero. 

^cyl 

RESULTS 

Figures 8 and 9 present a comparison of grid point displacement and element 
stresses from examples 1 and 2. The results are seen to be identica 1 from both 
these examples. The CPU time, however, on IBM 370/3031 was 750 seconds for 
example 1 as compared to 185 seconds for example 2. The expected behavior about 
the k = 2 natural frequency of 1814 Hz can be seen in both these figures. 

Figures 10 (k = 0), 11 (k = lc) and 12 (k = 2c) are from example 3. The 
k = 0 excitation consists of axial base acceleration and directly applied loads. 
The selected frequency band of excitation, 1700-1920 Hz, lies between the second 
out-of-plane disc bending mode frequency (1577 Hz, k = 0, Table 2) and the first 
in-plane shear mode frequency (1994 Hz, k = 0, Table 2). Since the excitation 
is parallel to the axis of rotation, only the former mode responds. The k = 1 
excitation is due to lateral base acceleration only. Although the frequency 
band of input base acceleration is 1700-1920 Hz, the rotation of the bladed disc 
at 600 Hz splits the input bandwidth into two effective bandwidths: 

(1700 - 600) = jnOO to (1920 - 600) = 1320 Hz, and 

(1700 + 600) = 2300 to (1920 + 600) = 2520 Hz. 

The only k = 1 mode in these effective bandwidths is the first torsional 

mode of the blade with the disc practically stationary (2460 Hz, k = 1, Table 2). 

This is shown by the out-of-plane displacement magnitudes of grid points 18 
(blade) and 8 (disc) respectively (Figure 11). The k = 2 excitation consists of 
directly applied k - 2c loads. The out-of-plane displacement magnitude of grid 
point 18 (Figure 12) compares well with that obtained in example 2 (Figure 8). 
Table 3 lists the out-of-plane displacement response of grid point 18 as obtained 
in examples 2 and 3. The marginal difference in response in example 3 is due to 
the Coriolis and centripetal acceleration effects at a rotational speed of 600 
revolutions per second. 

Results from examples 4 and 5 compare well with ti> from example 3, and 
are shown in Table 4. 


5. CONCLUSIONS 


1. Theoretical aspects of a new capability, developed and added to the 
general purpose finite element program NASTRAN Level 17.7 to conduct forced 
vibration analysis of tuned cyclic structures rotating about their axis of 
symmetry, have been presented. 

c> 

2. The effects of Coriolis and centripetal accelerations together with 
those due to the translational acceleration of the axis of rotation have been 
included. 

3. A variety of user options is provided to specify the loads on the 
rotating structure. 

4. Five interrelated examples are presented to illustrate the various 
features of this development. 
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APPENDIX A 

TERMS APPEARING IN EQUATION 4 


r u > 


{u} 


\->\ 
1 U 2 J 


(Al) 


r P. 


{P} * 


i 


(A2) 


{ V = 



(A3) 


global rn 1 
[MJ = [M 2 = [T GB J 


global 
[B] = [B 2 


m 

0 

0 


0 

m 

0 


0 

0 

m 


[T gb J 


(A4) 


(A5) 


global 


[B,] = [B, ] = [T GB J 


0 0 0 
0 0 -m 

0 m 0 


[T GB J 


(A6) 


global 

[Kj = [K 2 (includes all differential stiffness contributions) (A7) 


global Gb< T 


CM } 2 « CM, 2 * [T J 


0 0 
m 0 
0 m 


[T gb J 


(AS) 


Cm 2 1 * [t gb 3 


m 0 0 

0 me ms 
0 -ms me 


(A9) 
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^r> 



c 5 cos n t , s s sin J2 1 , and 



(A10) 

(All) 


s 
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APPENDIX B 

INERTIAL LOADS DUE TO BASE ACCELERATION 


The acceleration of the axis of rotation generates inertial loads at all 
grid points of the complete structure. In this appendix, the generation of these 
inertial loads and their transformation to frequency-dependent circumferential 
harmonic components are discussed. 

As given by equation (4), the inertial forces on the 
three translational degrees of freedom at an arbitrary point P of the modelled 
cyclic sector, expressed in the global (displacement) coordinate system, are 


where 


<P G ) = -[m 2 ]{R 0 > ■ [T BG ]{P B } , 



(Bl) 


(B2) 


with c = cos fit and s = sin fit. 

Since all the cyclic sectors are identical in all respects except for the 
specified loads, no generality is lost in assuming, for simplicity, that the 
modelled sector is the n = 1 sector. Equation (Bl) can, then, be rewritten as 


where 


<p g > ■ [t bg J 


1 

0 

0 



{P B } 



c 


n J 


, 


c R = cos (n-1 • 1 • 2tt/N), and 
s n = sin (rvT • 1 • 2tt/N ) 



Substituting equation (B3)in equations (13), and noting that 


(B3) 


(B4) 



= 0 


(B5) 
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or poe;; t. 


I* n = 0 


l c n • cos (fPT. k • 2tt/N) * N/2 , k s l 


0 , Ml , 


Js n • cos(rvT • k • 27t/N) = 0 , 


Jc • sinfrFT • k • 2tt/N) = 0 , 


Is n • sinCnHT • k • 2n/N) = N/2 , k=l 


= 0 


. kf 1 , 


(B5) 

(contd) 


the circumferential harmonic components of the base acceleration loads become 


tV B 


{P 0 } 


rnl C 


[I BG ] 0 


r 0 


B 


{p .c } . [ T DU j p 



( M k" = o) 


("k" = 1c) 


("k" = Is) , and 


{ pkc,ks } = {Q} t all other n k „ t 


(B6) 


*# M 


In the present development, the components of base acceleration X Q , Yq and 
Zq are considered to be sinusoidal of frequency o>, and are specified as 


x 0 ■ x 0, mag cos(u,t + V 


\ 


V 0 = v 0, mag cos <“ t + V • and 


2 0 ■ Z 0, mag cos(u,t + V 




(B7) 


J 
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From equation (B2) therefore, we can write 
P X B ‘ ~ m *0,mag. cos ^ t + * 

Py B * " m ^o,mag cos nt * cos ( wt + $y) + ^0,mag s1n * cos ( ut + and 
P Z B = -m[-Yo <ma g sinftt • cos(wt + <J> y ) + Z 0(fr)ag cos ft, cos(u)t + «|> Z )J. 

The cosine and sine products in equations (B8) can be expressed in terms of 
individual cosine and sine terms with frequencies (w + ft) and (to - ft), 

The following conclusions about base acceleration loads can, therefore, be 
drawn by substituting equations (B8) into equations (B6): 

•* n 

1. The axial component of base acceleration, Xq(w), contributes to P at exci- 
tation frequencies to. 

• • M 

2. The lateral components of base acceleration, Yg(to) and Z Q (to), contribute to 
F* c and F^ s at excitation frequencies (to ± ft) fur each to specified, 
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NOMENCLATURE 


B 

B V 

D 

G 


I» J. K 

A A /V 

V J B’ K B 

AAA 

i. j, k 

K 

k 


M 


m 

N 

P 


Q 

it 

R o 

r, R, p 
T 


t 

U 


u 

W 

8 


w 


Damping matrix 

Coriolis acceleration coefficient matrix 
Rayleigh's dissipation function 
"Symmetric Components" transformation matrix 
Unit vectors along Inertial XYZ axes 
Unit vectors -long Basic X B Yg Z B axes (Figure 1 ) 


Unit vectors along Global xyz axes J 
Stiffness matrix 
Circumferential harmonic index 
Time harmonic index 

Mass matrix, number of time intervals per period (Figure 3) 
Centripetal acceleration coefficient matrix 
Base acceleration coefficient matrix 
Mass 

Nun-bcr of cyclic sectors in the complete structure 
Load vector 

Aerodynamic coefficient matrix 
Base acceleration vector 
Position vectors (Figure 1) 

Kinetic energy, coordinate system transformation matrix 
Time 

Strain energy 

Physical displacement degrees of freedom 
Virtual work 
Rotational frequency 
Forcing frequency 
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NOMENCLATURE (Continued) 


Superscripts 


B • 

G 

K 

m 

n 

-0 

-Ac 

-As 

-kc ^ 
-ks 
-M/2 
-N/2 


Basic 

Global 

Independent solution set In "symmetric components" 
mth time instant 
nth cyclic sector 


Fourier coefficients ("symmetric components") 


NOTE : The above superscripts have been used in the paper to qualify various 

scalars, vectors and matrices, such as T, U, u. Pi M, K. They 
should not be confused with exponents. 
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ORIGINAL 



fljurt 3: Olrtctly Applied Ftrledlc loidt Sptelflod it 
functloni of Hint 


F1CUAE 3 : Orprill floxrtprt of Feretd *kr«tlpr- J»iiytlt pf Aem'rf 
Cyclic Structure! 































TAIlt I! MIKt’M. FttTWtS Of (lAFIPlt MdSICMS 


ORIGINAL PAC& U 

OF POOR QUALITY 


lit* 

ROttt lOfttl 

Rcc*l*r*t1vn 

Sp*fd 



A(F) - to»(K1T<2>2«/12) 

r, 

urTd pt.i 
Dor 

A{f) - tos|rCT*2-2»/12) 


B/2 C/ , 

A(f) 

A[t) • cot(p ; T , 2 , 2«/12) 


>«' z cyl. 

Alt) 




table ti irrtcT or corious ano centripetal accclzmtioss on Tn 
O lSHACCNtM RlS*0NSt or (RIO POINT It At too RPS, 


Fr«u»ncy 


1700 


1750 


1776 


1766 


1114 


1637 

1660 

1660 

1670 


7.2666 (-6/346.4 

7.6132 [-6/364.3 

1,3071 [-4/343,1 

! 1,3644 [.4/347,3 


2.1110 1-4/337.7 


3,4116 (.4/314,6 


4,1374 1-4/260, 6 


3.414} (-4/234,6 
2,1461 r-4/206,6 
1,2433 [-4/166,6 
7,6126 [-6/160.4 


2.3262 [-4/336.6 


3.7252 (.4/316,2 


4.617? [-4/266,6 


3.2(66 (-4/226.6 
2.0742 [-4/205. 3 
1.22H [-4/166.2 
7.535? [-6/564.3 



TAILE 2l ILAOED-DISC NATURAL FREQUENCIES 
Fr*ou*ncy (Moa* No.), N:, 

T ; .. — . | Nod* Dmrtptlon 


k • 0 4*1 

k • 2 

2)‘ ID 

208 (1) 

242 (1) 

591 (2) 

694 (2) 

622 (2) 

1577 (3) 

1123 (3) 

1B14 (3) 

2468 (5)** 

2460 (4) 

2433 (4) 



•Ml tn* Clreumftr«nt1i1 Mrmonlc Indl-x 

•• Nod* No. 4 for k * O *t 1664 H? r»pr*ttntt In In-pltnt tnm mod* nr-; *»c1t*d 
by tn* ipplltd forcti. 


tA6Lt 4 ; CONPARISON Of RESPONSE AT 1814 Hj, 


Srld Pt.Diip. or 
El*'. Strutts 


B «W 


16 (2 cy , t > 


11 • r «,1 


11 . 4 yy ,l 


11 • V.i 



*nt*ut»n 

Flour* 12s 2 ru . D1tp1*c*mtntl from Enimpl* 3, k • 2c, P1r»ctly Applltd lo*ds 


" * V.2 


* Tep tnd Bottom Fibres. 


In /Pnit*(d*sl 


6.4267 [-4/62,6 


4.9177 E-4/266.E 


1.4B41 [ 3/64,7 


2.0651 ( 2/63.4 


1,0774 [ 2/64,7 


1.4677 [ 3/263,3 


2.24E6 E 2/260,3 


1.6510 E 2/253,0 


5,4265 (-4/62. 6 


4,6160 [-4/2(6. 6 


/ [ 3/84,7 


2.0852 E 2/83.4 


1,0775 [ 2/64,7 


1.4676 E3/263.3 


2,2491 E 2/260.4 


1.6511 ( 2/253.0 


5 4255 E*4‘£2,6 


4.518C [.4/266,6 


1,4642 tJ’Ei.7 


1.0776 E2 '64.7 


1.4(76 [3/ 763.3 


2.24 


,6552 £2 253,0 


» 

















































































